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In paper I of this series we discuss how magnification bias distorts the 3D correlation function by 
enhancing the observed correlation in the line-of-sight (LOS) orientation, especially on large scales. 
This lensing anisotropy is distinctive, making it possible to separately measure the galaxy-galaxy, 
galaxy-magnification and magnification-magnification correlations. Here we extend the discussion 
to the power spectrum and also to redshift space. In real space, pairs oriented close to the LOS 
direction are not protected against nonlinearity even if the pair separation is large; this is because 
nonlinear fluctuations can enter through gravitational lensing at a small transverse separation (or i.e. 
impact parameter). The situation in Fourier space is different: by focusing on a small wavenumber 
k, as is usually done, linearity is guaranteed because both the LOS and transverse wavenumbers 
must be small. This is why magnification distortion of the galaxy correlation appears less severe 
in Fourier space. Nonetheless, the effect is non- negligible, especially for the transverse Fourier 
modes, and should be taken into account in interpreting precision measurements of the galaxy 
power spectrum, for instance those that focus on the baryon oscillations. The lensing induced 
anisotropy of the power spectrum has a shape that is distinct from the more well known redshift 
space anisotropics due to peculiar motions and the Alcock-Paczynski effect. The lensing anisotropy 
is highly localized in Fourier space while redshift space distortions are more spread out. This 
means that one could separate the magnification bias component in real observations, implying that 
potentially it is possible to perform a gravitational lensing measurement without measuring galaxy 
shapes. 

PACS numbers: 98.80.-k; 98.80.Es; 98.65.Dx; 95.35.+d 



I. INTRODUCTION 

The effect of magnification bias on the 3D 
galaxy/quasar correlation function was studied in 
paper I [![. (Galaxy and quasar can be considered 
synonymous hereafter.) With the important exception 
of the classic paper by Matsubara , previous work on 
how magnification bias modifies clustering observations 
has largely focused on the 2D angular correlation 
function 3, 1 i, H, 0, i, i, M MMm- The novelty 
of the 3D correlation function, as emphasized by 0] and 
is that magnification bias makes it anisotropic. In 
this paper, we extend our previous analysis by studying 
the anisotropy in Fourier and redshift space. At first 
sight, the extension to Fourier space might seem a trivial 
exercise. The calculations are indeed straightforward, 
but as we will see, the results are far from obvious: there 
are important qualitative differences between the results 
in Fourier space and real space that go beyond the usual 
wavenumber-position (k-x) duality. 

Let us recall the situation in real space, as depicted 
in paper I. The anisotropy of the observed 3D correla- 
tion function can be understood intuitively as follows. 
The correlation function is measured by pair counts of 
galaxies. A pair of galaxies that are aligned along the 
line-of-sight (LOS) behave differently from a pair ori- 



ented transverse to the LOS. In the former case, the 
closer galaxy can lens the background one. The same 
does not happen in the transverse orientation. The net 
effect is an anisotropy in the observed correlation func- 
tion, induced by gravitational lensing (or equivalently, 
magnification bias; we refer to this effect as magnifica- 
tion distortion). This reasoning suggests the gravita- 
tional lensing corrections are largest for a pair of galax- 
ies oriented along the LOS. Indeed, the corrections can 
be quite significant: consider for instance a LOS sepa- 
ration of ~ 100 Mpc/h; the intrinsic galaxy correlation 
is rather weak on such a large scale, but the lensing in- 
duced correction can be quite substantial, since for the 
LOS orientation, the relevant lensing impact parameter, 
i.e. transverse separation, is zero (keep in mind also that 
the lensing effect grows with the LOS separation while 
the intrinsic galaxy correlation generally drops with sep- 
aration). In other words, for the LOS orientation, scales 
that otherwise would be considered linear can in fact be 
secretly affected by nonlinear fluctuations via lensing. A 
large separation |<5x| = 5x^ + |<5x^p does not guaran- 
tee linearity because nonlinear fluctuations can sneak in 
through lensing with a small transverse separation |(5x^ |. 

This peculiar mixing of linear intrinsic galaxy fluctua- 
tions with nonlinear lensing fluctuations does not arise in 
Fourier space. A small net wavenumber |k| is sufficient 
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to guarantee that both the LOS component fc|| and the 
transverse component |kj^| are smah. Nonhnear lensing 
corrections cannot sneak in as long as one focuses on a 
small |k|, as is usually done. This immediately tells us 
that the anisotropy of the observed galaxy correlation 
must appear milder in Fourier space. Our primary goal 
here is to quantify this. 

As discussed in paper I, the magnification bias induced 
anisotropy in the observed 3D galaxy correlation has two 
implications. First, precision measurements of the galaxy 
correlation must take into account such magnification 
distortion. These include future galaxy surveys that hope 
to determine the baryon oscillation scale to high accu- 
racy [11, [H m, [O, 111 [H iOl- Second, the distinctive 
anisotropy pattern makes it possible in principle to sep- 
arately measure the galaxy-galaxy, galaxy-magnification 
and magnification-magnification correlations. (The last 
correlation was generally ignored in previous papers that 
focused on angular correlation between galaxies at widely 
separated redshifts [1, H, [H [HI, [H, [II] , where the galaxy- 
magnification correlation dominates.) Achieving such a 
separation requires that one understands other sources 
of anisotropy. We therefore extend our Fourier analysis 
here to incorporate the anisotropy due to both peculiar 
motions and the Alcock-Paczynski effect. 

The rest of the paper is organized as follows. In SjlTl we 
derive and numerically compute the magnification distor- 
tion of the observed galaxy clustering - ^11 Al summarizes 
the results from paper I on the correlation function while 
miB\ focuses on the power spectrum. Redshift space dis- 
tortion due to peculiar motion is next incorporated in 
mm and we conclude in mVl In Appendix [XI we discuss 
the Alcock-Paczynski anisotropy. 

Before we start, it is useful to point out several re- 
lated papers. Vallinoto et al. [2l| explored the impact 
of lensing, especially magnification bias, on the baryon 
oscillation signal in the real space correlation function. 
Their results are consistent with ours in paper I, though 
they focus exclusively on pair separations that are ori- 
ented transverse to the LOS, and their work is therefore 
more connected to our paper on the angular correlation 
function [13] . Wagner et al. [1^ examined the anisotropy 
of the 3D correlation that is introduced by light cone ef- 
fects. A discussion of the classic paper by Matsubara [§] 
can be found in paper I. Both |2] and paper I focused on 
the real/configuration space correlation function, though 
peculiar motions and the Alcock-Paczynski effect are also 
treated in ^]. A recent paper by Zhang & Chen j2^] ex- 
plored the effects of gravitational lensing in Fourier space 
in the context of supernova observations. LoVerde et al. 
[25| examined the impact of magnification bias on inte- 
grated Sachs- Wolfe measurements. 



overdensity Jobs: 



Sohs = + 6n 



(1) 



which is a function of the galaxy position, specified for 
instance by the radial comoving distance x and the angu- 
lar position 6. The magnification bias correction is given 
by Biii[3,i: 

Sf, = (5s - 2)k (2) 
where k is the lensing convergence: 



k{x, e)= dx 
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(f> is the gravitational potential, and is the 2D Lapla- 
cian in the transverse directions. We assume a flat uni- 
verse - generalization to an open or a closed universe is 
straightforward. The symbol s stands for 



dlog^oN{< m) 
dm 



(4) 



where iV(< to) is the cumulative number counts for 
galaxies brighter than magnitude to. This assumes the 
galaxy sample is defined by a sharp faint-end cut-off. A 
broader definition of s for a more general galaxy selection 
is given in Appendix A of paper I. 

We define the galaxy bias b by dg — bS, where S is the 
mass overdensity. Eq. ([!]) can then be rewritten as 

Jobs c , 5s - 2 



(5) 



The relative importance of the intrinsic clustering and 
the magnification bias correction is therefore controlled 
by, among other things, the sample dependent ratio (5s — 
2)/&. 

The precise values of s and b depend sensitively on 
details of how the galaxy/quasar sample is selected, for 
instance subject to color cuts and so on. Unless otherwise 
stated, we adopt throughout this paper the value (5s — 
2)/6 = 1 to illustrate the effect of magnification bias on 
clustering measurements (see paper I for more details). 
An implicit assumption is the linearity of the galaxy bias, 
a subject we will return to in mV\ 

In all illustrative examples below, we employ the fol- 
lowing cosmological parameters: the Hubble constant 
h = 0.7, matter density ri„i — 0.27, cosmological con- 
stant fl\ — 0.73, baryon density fib — 0.046, power spec- 
trum slope n = 0.95 and normalization erg = 0.8. We 
employ the transfer function of [2^ , and the prescription 
of [23] for the nonlinear power spectrum. In all equations 
we use units where the speed of light is unity: c = 1. 



The Correlation Function 



II. MAGNIFICATION DISTORTION 

Given an intrinsic galaxy overdensity dg, magnification 
bias introduces a correction (5^ to the observed galaxy 



Here, we summarize the main results of paper I. Includ- 
ing lensing magnification, the observed two-point corre- 
lation function is given by: 

Cobs(Xl,^i;X2,02) = (Jobs(Xl,^l)'^obs(X2,02)) (6) 
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= CsslXl, ^i; X2, ^2) + X2, ^2) 

+£.gp.{X2, ^2; Xl, ^1) + ^A'A'(Xl, X2, ^2) 

where the magnification bias corrections, the galaxy- 
magnification and magnification-magnification correla- 
tions, are: 

?gM(Xl,^i;X2,^2) +CgM(X2,02;Xl,^l) = (7) 

^i?2f]„(5s-2)(l + z)|x2-Xi| 



Cmm(xi, ^i; X2, 02) - [-Hl^lr,,{bs - 2)f (8) 



dx' 



(x - x')x' 



X 



(1 + ^0 



/\2 



Here we have Taylor expanded xi ^-nd X2 around the 
mean x retained the lowest order contributions. The 
intrinsic (unlensed) galaxy auto-correlation, or galaxy- 
galaxy correlation, is 



Css(xi,0i;x2,02) = 
i,aW{xi^X2? + noi-e2f) - 



(9) 



{2n) 

ignoring for now the issue of redshift distortion, which 
will be addressed in §1111 Note that xi and X2 refer to 
the points corresponding to xi, 61 and X2, ^2- Note also 
Pmm, Pgm and Pgg dcuotc respectively the mass- mass, 
galaxy-mass and galaxy-galaxy power spectra. 

The observed correlation function is a sum of all three 
correlations above (eq. [7], [5] and [5]). (A discussion 
of their higher order Taylor corrections can be found 
in Appendix B of paper I.). Viewed in this way, the 
anisotropy of the lensing induced corrections is quite 
striking: 2) -|-^g^(2, 1) scales linearly with the line- 

of-sight (LOS) separation |x2— Xi| (i-^- it increases rather 
than decreases with the separation!), and is indepen- 
dent of the LOS separation. The intrinsic galaxy auto- 
correlation is isotropic and generally decreases with 
separation. 

We can summarize the distinctive lensing induced 
anisotropy in the observed correlation function as follows: 



£.ohs{5x^5x±_) 



S.gg{^ 5x^ + 5x1) (10) 
+f{5x\_)5x + g{5x\_) 



where 5x and 5x±^ are the LOS and transverse 
separations respectively, f5x represents the galaxy- 
magnification correlation and g represents the 
magnification-magnification correlation. Here, / 



and g are functions of the transverse separation only, 
and are determined by the galaxy-mass and mass-mass 
power spectra. This distinctive form of the anisotropy 
allows us in principle to separately measure ^gg, f and g, 
from which we can infer the galaxy-galaxy, galaxy-mass 
and mass-mass power spectra. For instance, at any given 
Sxj_, plotting ^obs as a function of the LOS separation 
Sx would reveal a linear contribution at sufficiently large 
(5x's where £,gg is very small. Its slope tells us / and its 
extrapolation to 5x = tells us g. Subtracting Sxf + g 
from ^obs then yields ^gg. This is illustrated in Fig. 2 
of paper I, where we also present order of magnitude 
estimates for the ratios ^^^/^gg and Cg^/^ss' 
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(11) 



The symbol A^(/c) denotes the dimensionless variance at 
scale k and redshift z: ^tt Pmm{k) / {2t:Y . Here, fc** ~ 
1 / -^/^xM-^a^, while fc* is equal to either l/5x^^ or fcm, 
whichever is smaller {km is the scale where k^Pmm{k) 
peaks; fcm > 3 h/Mpc). These estimates work reasonably 
well except for separations around or beyond the zero- 
crossing scale. 

For the LOS orientation, where A^(fc,) can be much 
larger than A^(fc**), the magnification bias corrections 
and ^g^i can dominate over the intrinsic clustering 
correlation ^gg. The implications for baryon acoustic 
oscillation measurements are summarized at the end of 

am 



B. The Power Spectrum 

For surveys with a simple geometry, the power spec- 
trum is often the more popular quantity to measure. 
Suppose the survey geometry is specified by T/F(x) = 
M^ll (x||)Wj_(x_l), where x specifies a location with x\\ be- 
ing the LOS component and x^ the transverse compo- 
nent. For instance, a top-hat geometry in the radial di- 
rection is described by 

W\\{x\\) = 1//L if ~L/2<x\\<L/2 (12) 
~ otherwise 

with L being the radial span. A Gaussian geometry in 
the radial direction would be described by 

W^||(a;||) = (7r(T2)-i/4exp[-a;f/(2a2)] (13) 

Our normalization convention is that J d^xW^ = 1. Note 
that the origin = is chosen to be located at the mean 
redshift of interest. 

The Fourier counterparts of eq. (O and ([5]), taking 
into account the effects of the window, are 



2Pg^(k) = -Hl^lm{5s -2){l + z) 



(14) 
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G(fc||) 



and 
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P^.ik) = [-i?o^f]™(5s-2)]^|W-||(fc||)|^ (15) 



(27r) 



where W\\ and W_l are Fourier transforms of the windows 
and Wj_, and G is defined as foUows: 

G(fc||) = /da;ida;2|a;i-:E2|e*'="("^-"^)W^||(xi)VK||(x2ll6) 



where xi and X2 represent the LOS distance. 

The galaxy power spectrum is windowed in the usual 
way: 

(Pk' 



(27r)3 99 



Pll-%k')\W{k-k!r (17) 



where W is the Fourier transform of the total window W ^ 
and ^^3™° represents the true/unwindowed galaxy power 
spectrum. 

The observed power spectrum is the sum: 

Pobs(k) = Pgg{k) + 2Pg,,(k) + P^^(k) (18) 

where Pobs(k) is defined to be 

Pobs(k) = j d3a;id3x2(5obs(xi)(5obs(x2)) (19) 
W(xi)T/F(x2)e^''-(''i-^^) 



The survey window enters into these three contribu- 
tions to the observed power in distinct ways. For the 
galaxy-galaxy power spectrum, the window function is 
convolved as usual under an integral with the true 3D 
power. For P^^i, only the transverse window function 
is convolved under an integral with the power spectrum 
{Pmm)- The LOS window function is not convolved un- 
der an integral at all: P^^ is directly proportional to 
|W^llp. For Pg^, it is also directly proportional to some 
generalized LOS window function G(fc|| ). This interesting 
behavior of the magnification-magnification and galaxy- 
magnification power spectra can be traced to the unique 
anisotropics of their real space counterparts: referring 
back to eq. PU)) . the appearance of \W\\\'^ in P^^ is re- 
lated to the fact that g is independent of the LOS sepa- 
ration, and the appearance of G in Pg^ is related to the 
linear dependence of f6x on the LOS separation. 

The precise form of the window functions depends on 
the exact geometry. For the example of the top-hat win- 
dow (eq. [H]), they are 




10-3 10-2 
k„[h/Mpc] 



FIG. 1: The multiplicative LOS windows \W\\f and Gy as a 
function of the LOS wavenumber fcy . The dashed lines are 
for a top-hat geometry, with L — 706.6 Mpc/h and the solid 
lines are for a Gaussian geometry, with a = 204 Mpc/h. 



sin (kuL) — cos (fciiL) — 1 

/c||L " " 



At low ~ L while G - 

The corresponding expressions for the example of the 
Gaussian window (eq. [13| ) are 



|W^||(fc||)l' = V4 
G(fc||) =4a2 



1 -cr^fe; 



(21) 



dye 



At low fc|| , I p ~ y/Ina, while G ~ 4cr^. Illustrations of 
1 14^11 p and G can be found in Fig. [l]for both the top- hat 
and Gaussian geometries. 

For completeness, let us also give the window func- 
tions in the transverse directions for a top-hat (a circle 
of radius R) : 



W^iki_) 



l/\ArP2 if x±<R 
otherwise 

drrJo (r) 



(22) 
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klR\ 


Jo 



k± 



Ji(fc±P) 



where Jo and Ji are the Bessel functions, and for a Gaus- 



(20) 



W^{x^) - (^a2)-i/2exp [-4/(2^2)] 



(23) 
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(a) (b) 

FIG. 2: z = 1: (a) Contours of constant P^hs (solid) and Pgg (dotted), left to right: logio (P/[Mpc/li]^) = 4.02 (red; double 
contours), 3.93 (blue: double contours), 3.83, 3.73, ... 3.33 (black), (b) Various power spectra normalized by the same BBKS 
(no baryon) galaxy power spectrum: Pobs for fcy = (red solid), Pgg (black dashed), monopole of Pobs ((5s — 2)/6 = 1 for red 
dotted and (5s ~ 2)/b — 2 for red dot-dashed). Note that fc^ = fcy + fc^. A Gaussian window is assumed with a = 204 Mpc/h 
for both panels, and (5s — 2)/6 = 1 is adopted throughout except for the red dot-dashed curve. 



W±{k±) = VAttct^ exp 

Note that in Fourier space the galaxy power spectrum 
Pgg can itself be anisotropic if the window function is 
anisotropic. In the Gaussian case, if the same cr were 
chosen for both W± and W\\ , the windowed galaxy power 
spectrum will remain isotropic. In our computations be- 
low in this section, we adopt this special choice in order to 
more clearly show the anisotropy induced by magnifica- 
tion bias. We have checked that using the top- hat geom- 
etry yields a rather similar lensing anisotropy as long as 
one makes the choice 2R ~ L ~ ^/l2a (the latter equal- 
ity is chosen such that the Gaussian G and the top-hat 
G have the same low fcy limit). 

As in the case of the correlation functions, it is useful to 
give order of magnitude estimates for the ratios Pfi^/ Pgg 

and Pgf,/Pgg: 
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(24) 



where fc^ = fc^ -I- fc^ . The above expressions are approxi- 
mate. For instance, we have approximated Pmm(k±x/x') 
(eq. [13]) by Pmm{k±)- Nonetheless, they agree with the 
exact numerical integration to within factor of a few, and 
they illustrate several important points. 

The presence of and G means that the effects 



of magnification bias are largest for low k\\ 's. The sim- 
plest limit to consider is the one with k± ^ fc|| (i.e. a k 
vector that is oriented transverse to the LOS) in which 
case the factors of Pmm cancel out in the ratios, since 
k ^ k±. In the smah /c|| limit, |I^||P ~ L and G ~ 
where L is the width of the redshift bin over which one is 
measuring the power spectrum. One can see that for in- 
stance at z ~ 1.5 where xHq ~ 1, and for (5s ~2)/b ^ 1, 
we have P^^/P^g - O.lHoL and 2Pg^,/Pgg - 0.4(i7oi)^ 
A choice of L ~ 900 Mpc/h (which corresponds to the 
redshift interval 1.5 ± 0.35) yields P^^/Pgg ^ 0.03 and 
"^Pgfi/Pgg ~ 0.04. lu othcr words, the total effect of 
magnification bias is quite modest, ~ 7% in this config- 
uration. It should be kept in mind that (1) this estimate 
increases with L, and (2) the P^^/Pgg ratio increases 
strongly with redshift: the cubic dependence on {xHq) 
and quadratic dependence on 1 + z means this ratio rises 
rapidly beyond z ~ 1.5. 

Nonetheless, it is perhaps a little surprising that mag- 
nification bias appears to have a much more modest effect 
on the observed clustering in Fourier space compared to 
real space, where the corresponding ratios ^^^/S,gg and 
2S,gfj,/^gg cau Tcach order unity or even higher in the LOS 
orientation (eq. [11]; see paper I for details). The funda- 
mental reason is the absence in Fourier space of this po- 
tentially large boost factor, A^(A:*)/A^(A;*») in eq. (fTTj) . 
that is present for the correlation function. Consider 
a separation vector (5x that is oriented along the LOS, 
which is the orientation that maximizes the magnification 
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(a) (b) 

FIG. 3: Analog of Fig. [2]for z = 1.5. The contours in (a) are, left to right: logio (P/[Mpc/h]^) = 3.86 (red; double contours), 
3.76 (blue: double contours), 3.66, 3.56, ... 3.16 (black). A Gaussian window is assumed with a = 270 Mpc/h. 




(a) (b) 

FIG. 4: Analog of Fig. [2]for z = 2. The contours in (a) are, left to right: logio (P/ [Mpc/h] ^) = 3.71 (red; double contours), 
3.61 (blue: double contours), 3.51, 3.41, ... 3.01 (black). A Gaussian window is assumed with a = 323 Mpc/h. 



bias effect. If \6k\ is sufficiently large, the galaxy-galaxy 
correlation ^gg is essentially determined by the power 
spectrum in the linear regime (low fc*, ^ l/l'^xj in eq. 
|llj ) and is quite weak, while the galaxy- magnification 
correlation and the magnification-magnification cor- 
relation are sensitive to the power in the nonlinear 
regime (high ^ l/6x± in eq. [11]) and can be appre- 



ciable. In other words, in the LOS orientation in real 
space, one is comparing intrinsic galaxy fluctuations and 
lensing fluctuations on very different scales. 

Consider, on the other hand, a vector k that points 
in the transverse direction, which is the Fourier analog 
of a LOS (5x. In this case, k ^ k± and factors of the 
mass power spectrum simply cancel out of the ratios in 
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(b) 



FIG. 5: Analog of Fig. [2]for z = 3. The contours in (a) are, left to right: logio (P/[Mpc/h]^) = 3.47 (red; double contours), 
3.37 (blue; double contours), 3.27, 3.17, ... 2.77 (black). A Gaussian window is assumed with g — 397 Mpc/h. 



eq. p4|) . There is no boost coming from a ratio of powers 
on very different scales, as in the case of the correlation 
function. In other words, the mixing of magnification 
bias corrections with the intrinsic galaxy clustering term 
occurs in a very different manner in Fourier space than 
in real space. 

Ultimately, the correlation function and the power 
spectrum are related by Fourier transform, and so should 
really contain the same information. The important 
point to keep in mind, however, is that for cosmological 
purposes, one often focuses on scales that are perceived 
to be linear: in the case of the correlation function, that 
means a large |^x|, and in the case of the power spec- 
trum, that means a small |k|. With the presence of mag- 
nification bias, a large |(5x| no longer protects one against 
nonlinearity i.e. a separation (5x pointing along the LOS 
(i.e. small |i5x_l|) is subject to large magnification correc- 
tions even if |(5x| is large. On the other hand, a small |k| 
means both and |k^| must be small, which protects 
one against nonlinear fluctuations. 

Let us return to the order of magnitude estimates in 
eq. I|24p. to see what happens if one considers angular 
averages. The monopole of any power spectrum P is 
defined as 



monopole of P(k) 



tt/2 



P(k) sinOkdek 



(25) 



where 6k is the angle between the LOS and k. Let us 
focus on the monopole Pfifi/Pgg for a Gaussian window, 
where the calculation is the simplest. Suppose one is 
interested in a scale k such that ka > 3. Because of 

the Gaussian in iWip, the integral over angle will be 



dominated by 9k ^ 7''/2, and so one can approximate 
Pmmik±) by Pmm{k). The remaining integral over angle 
is simple to do, and yields the ratio: 



monopole of 
monopole of P, 



gg 



5s 



which is valid only for /ccr > 3. Note how the width of 
the window a completely disappears from this ratio. A 
similar expression holds for the tophat case as well. This 
ratio is small for k/Ho ^ 1, unless one goes to a suffi- 
ciently high redshift. The corresponding ratio for galaxy- 
magnification cannot be worked out analytically because 
the window G has a more complicated form. On dimen- 
sional grounds, one expects this ratio to scale with a. In 
practice, we find that for redshifts where magnification 
bias matters, the monopole of the galaxy-magnification 
power spectrum is quite a bit smaller than that of the 
magnification-magnification power spectrum, in part be- 
cause of cancellations that occur under the angular aver- 
age. 

The intuition gained above from the order of magni- 
tude estimates (eq. [24], [26]) is confirmed by the exact 

-f^M ^''^d Pgg according to 
This is shown in Fig. [U - [U 
for redshifts z = 1, 1.5, 2, 3. The magnification distortion 
increases with redshift and, as discussed before, has the 
most noticeable effects for small fcy 's. Exactly how small 
k\\ needs to be to see a substantial effect depends on the 
LOS width of the survey/sample selection function. For 
a Gaussian window function the region of large magnifi- 
cation distortion is fc|| < cr~^. Panels (b) of Fig. O-IH] 
show that the monopole of the power spectrum is visi- 



numerical evaluation of P, 
eq. ([H]), ^ and ^ 
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bly distorted as well, the effect at low fc's is particularly 
severe at redshifts > 2. 

There is the interesting question of how magnification 
bias impacts baryon oscillation measurements. This was 
addressed in paper I for the real space correlation func- 
tion. Briefly summarizing: we found that the observed 
baryon acoustic scale can be shifted by up to ~ 3% in 
the LOS orientation, and up to ~ 0.6% in the monopole, 
depending on the exact values of the galaxy bias, red- 
shift and number count slope. The corresponding shifts in 
the inferred Hubble parameter and angular diameter dis- 
tance, if ignored, could significantly bias measurements 
of the dark energy equation of state (by up to ^ 15%). In 
Fourier space, there are several wiggles, and the magnifi- 
cation bias induced shift in the baryon oscillation scale is 
likely more sensitive to exactly how this scale is extracted 
from data. We therefore do not attempt to investigate 
this further in this paper. Given the earlier discussions, 
it is reasonable to expect that baryon acoustic oscillation 
measurements are less affected by magnification bias in 
Fourier space. However, it is worth emphasizing that 
magnification bias introduces scale and orientation de- 
pendent corrections to the observed power spectrum, and 
these corrections depend on uncertain factors such as the 
galaxy bias. The question is whether, in fitting the ob- 
served data for the baryon oscillation scale, one should 
introduce additional fit parameters to account for mag- 
nification bias, and what impact they might have on the 
measurement accuracy of the Hubble parameter and the 
angular diameter distance. This certainly deserves more 
study. 

It is also worth noting that, as can be seen from panels 
(b) of Fig. [2]-[5l the radiation-matter equality peak loca- 
tion/shape around k ^ 0.01 h/Mpc is likely significantly 
affected by magnification bias, and one must be careful 
in using it as a standard ruler [2^. 



III. INCORPORATING REDSHIFT 
DISTORTION DUE TO PECULIAR MOTIONS 

The observed redshift of a source galaxy is dependent 
upon both the distance to the source (the cosmological 
redshift) and the peculiar velocity of the source. In red- 
shift space, the observed galaxy density is (to first order 
in perturbations): 



where Sg and (5^ are as in eq. (H]), and 5y is 



Sv 



(1 + z) dv\\ 
H{z) dx\\ 



(27) 



(28) 



where U|| is the LOS peculiar velocity, and H{z) is the 
Hubble parameter at the mean redshift z. 

The observed two-point correlation function, instead 
of eq. ([ni, is now given by 

Us(l; 2) = igg{l; 2) + ig^{\- 2) + ig^{2- 1) (29) 



-HC^^(1; 2) + ^3.(1; 2) + iyg{2- 1) -I- C™(1; 2) 

where we have used the arguments 1 and 2 as the short- 
hand for the positions of the two points of interest in 
redshift-space. We have used the Limber approximation 
which makes the velocity-magnification cross-terms van- 
ish (the derivative with respect to a;|| in the velocity term 
pulls down a factor of fc||, which vanishes under the Lim- 
ber approximation; see [3]). The galaxy- velocity cross- 
correlation and the velocity auto-correlation are given by 
the well-known results of Kaiser [291 : 



^9.(1,2) = 
and 
U(l,2) = 



a D' 



a D' 



(27r)3 P^^™ 



(fc), 



3jk-(xi— X2) 



(30) 



d^k kl 
(27r)3 fc4 



k-(xi-X2) 



(31) 



where Xi — X2 has a LOS component Xi — X2, and trans- 
verse components of x(^i ~ ^2)- Here, a is the scale 
factor, D is the linear growth factor, and a' and D' are 
their derivatives with respect to conformal time. All time 
dependent quantities are evaluated at the mean redshift 
z. 

The observed power spectrum, instead of eq. (jlSp . is 
now given by 

Fobs (k) = Pgg (k) + 2Pg^ (k) -f P^^ (k) (32) 

+2Pg,(k) -f P„,(k) 
where the first line is as before (eq. [l3], [IS] & [HDj and 



Ps.(k) = 
and 

Pvv (k) — 



a D' 



a D' 



d^k' k\^ 



(27r)3 k 



d^k' k\i' 
(2^)3 k'' 



^P,„(fc')|M^(k-k')P(33) 



'Pmm 

(fc')r(k-k')|^34) 



To gain some intuition about the various effects at 
work, it is useful to adopt the following approximation: 
integrate out all the convolving windows as if they are 
delta functions. (This approximation is made only in 
this section, not in previous sections.) For instance: 



d^k' 
(2^ 



p(fc')|vK(k-k')r-p(fc) 



(35) 



where P represents Pgg, Pgm, Pgmk'^i^ /k' and so on, and 

the convolving window W could also be W± . This works 
well if the power spectrum is sufficiently smooth and the 
convolving window is sufficiently narrow. With this ap- 
proximation, we obtain 



Fobs (2, k) = Pgg{z,k) 



7 9 \ 2 

Id ^11 \ 



b k^ 



-Kcrt 



(36) 
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kJh/Mpc] k [h/Mpc] 



(a) (b) 

FIG. 6: Analog of Fig. |4]for z = 2, except redshift distortion is incorporated with fo/b = 0.475, and (Tz = 3 Mpc/h (eq. 
[36]). (a) Contours of constant Pobs with magnification bias (solid) and without magnification bias (dotted), left to right: 
logio (P/[Mpc/h]=*) = 3.71 (red; double contours), 3.61 (blue: double contours), 3.51, 3.41, 3.31, 3.21 (black), (b) Various 
power spectra normalized by the same BBKS (no baryon) galaxy monopole power spectrum: Pobs with magnification bias for 
fc|| = (red solid), Pobs without magnification bias for fcy = (black short-long dashed), monopole of Pobs without magnification 
bias (black uniform dashed), monopole of Pobs with magnification bias ((5s — 2)/ti = 1 for red dotted and (5s — 2)/b = 2 for red 
dot-dashed). Note that — k^^^-'r k\. A Gaussian window is assumed with a = 323 Mpc/h for both panels, and (5s — 2)/& = 1 
is adopted throughout except for the red dot-dashed curve. 



q{l + z)G{h) 



g2|#„(fc„)|2 



'dx'(x-x')^(i + ^')^- 



.{z',k^x/x') 



where 



d\nD 



(37) 



(38) 



evaluated at the mean redshift z, and 

and we have assumed a linear galaxy bias b (at redshift z). 
We have introduced an exponential factor exp[— ^jja^] 
which accounts for a possible (Gaussian) dispersion in 
redshifts, such as from photometric redshifts. Here Uz is 
the dispersion expressed in comoving Mpc/h (not in red- 
shift). One can also think of this exponential factor as 
modeling the effect of nonlinear or virialized peculiar mo- 
tions, though this description is at best approximate [30| . 
Strictly speaking, with a non-zero az, the multiplicative 
window G in eq. (j36p should be replaced by 

dzidz2dz[dz'2 , 



Jkn (Z1-Z2) 



27rcr? 



(39) 



W\i{zi)Wi\{z2) exp 







[ (4-^2)^1 


[ 2a2 J 


exp 


[ 2a2 J 



where zi, Z2, z'l, z'2 denotes LOS distances, not redshifts. 
We find that as long as the LOS width of the sur- 
vey/sample selection is large compared to az (i.e. 
L ^ az for a top-hat geometry, or a ^ az for a Gaus- 
sian geometry; see eq. [20] & [H]), the above Gz is well 
approximated by G as defined in eq. 

The first line of eq. excluding the expo- 

nential factor, represents the classic Kaiser distortion 
due to coherent infall. The second and third lines 
come from the galaxy-magnification and magnification- 
magnification correlations respectively. We illustrate all 
these effects in Fig. [B] and [7] by evaluating eq. for 
z = 2. Fig. [6] uses 0-2=3 Mpc/h, which corresponds to 
a velocity dispersion of ~ 300 km/s, or a redshift disper- 
sion of 0.003. This is the level of dispersion one expects 
from virialized motions on small scales. 

For the scales shown in Fig. [6l the Kaiser distortion 
dominates over the finger-of-god effect due to virialized 
motions: the contours of constant Pobs are elongated 
in the LOS direction for fe's larger than the radiation- 
matter equality scale, and compressed otherwise. It is 
also worth noting that magnification distortion survives 
the Kaiser distortion. The two kinds of distortions have 
fundamentally different shapes. Magnification distortion 
is localized to small fc|| 's, whereas the Kaiser distortion is 
more spread out. The localized nature of magnification 
distortion originates from the fact that the corrections 
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(a) (b) 

FIG. 7: Analog of Fig. [g] except that = 30 Mpc/h. The contours in (a) are, left to right: logio (P/[Mpc/h]^) = 3.71 (red; 
double contours), 3.61 (blue: double contours), 3.51, 3.41, ... 1.01 (black). In panel (b), note how the monopole (red dotted 
with magnification bias and black uniform dashed without) is higher than the fey = power (red solid with magnification bias 
and black short-long dashed without) for small fe's due to the Kaiser effect, and lower for high fc's due to the finger-of-god effect 
(from the large CTz). 



it introduces (the second and third terms on the right 
hand side of eq. |36| ) are proportional to the multiplica- 
tive windows G(fc||) and |W^||(fc||)P, both of which peak 
at small fc||'s. 

Interestingly, the Kaiser distortion vanishes at fc|| = 0, 
exactly where magnification bias has the largest effect. 
This suggests it should be possible to disentangle the 
two different distortions from data. For instance, the one 
free parameter that controls the Kaiser distortion, 
can be determined from the observed power spectrum 
anisotropy by excluding from consideration the small /e|| 
modes. With this in hand, one should be able to predict 
the dotted contours such as those in Fig. [5^. The dif- 
ference between the observed (solid) contours and the 
dotted ones then gives us the magnification bias cor- 
rections. The galaxy-magnification and magnification- 
magnification contributions can be further separated 
from each other by using the distinctive shapes of their 
respective multiphcative windows G and (for 
instance, the latter is positive definite whereas the former 
can go negative; see Fig. [T|). 

Fig. [7] is analogous to Fig. [H] except that Uz is in- 
creased to 30 Mpc/h, corresponding to a velocity disper- 
sion of ^ 3000 km/s, or a redshift dispersion of 0.03. 
This larger value for Gz is chosen to mimic the effect of 
photometric redshifts. In contrast to the case of cr^ = 3 
Mpc/h, one can clearly see here a finger-of-god effect of 
sorts: the contours of constant Pobs are compressed in 
the LOS direction, for sufficiently large fc's. As before, 



magnification distortion is clearly visible, being well lo- 
calized to small fc|| 's. Its distinctive shape makes it in 
principle distinguishable from both the Kaiser distortion 
and the finger-of-god effect, which are more spread out 
on the fcj^ — fc|| plane. 

It is also worth noting that in cases where the power 
spectrum is strongly anisotropic, such as in Fig. [71 the 
monopole is probably not the most relevant quantity to 
consider. This is because some orientations are much 
more noisy than the others, and one might not want 
to weigh them equally (or more precisely, according to 
eq. [23 )• In general, with sufRcient redshift accuracy, 
one should make use of the full 3D information avail- 
able. In this vein, it is not uncommon to consider higher 
multipoles of the anisotropic power spectrum [3l| . This 
is especially useful for analyzing the Kaiser distortion 
since it gives rise to only 2 extra multipoles. However, 
a multipole expansion is likely not helpful in analyzing 
magnification distortion, due to its localized nature. It 
is probably more useful to take advantage of the spe- 
cial dependence on G{k\\) and |W^||(fc||)P of the galaxy- 
magnification and magnification-magnification terms (eq. 
|36|). How to optimally extract these two contributions 
from noisy data deserves further study. 

In addition to redshift distortion due to peculiar mo- 
tions, another well known effect is the so called cosmolog- 
ical distortion, or Alcock-Paczynski effect [H, [s^. This 
is taken up in Appendix W\ 
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IV. DISCUSSION 

In paper I, we examined the effects of magnification 
distortion in real/configuration space, and here we have 
extended the analysis to Fourier and redshift space. The 
observed galaxy/quasar correlation is endowed with a 
distinctive lensing induced anisotropy. This is encap- 
sulated in eq. pU|) for real space and eq. (fH)) . 
P?)) and (HH]) for Fourier space: the linear dependence 
on the LOS separation Sx in the real space galaxy- 
magnification correlation gives rise to the multiplica- 
tive window G in the galaxy-magnification power spec- 
trum. Likewise, the independence of the real space 
magnification-magnification correlation on Sx accounts 
for the appearance of the multiplicative window p in 
the magnification-magnification power spectrum. 

Qualitatively, the galaxy correlation becomes en- 
hanced for the transverse modes in Fourier space, and 
for pairs oriented along the LOS in real space. Quantita- 
tively, the degree of enhancement is rather different for 
the dual spaces. As explained in and ijllBl magnifi- 
cation distortion is less severe in Fourier space: as long 
as one focuses on modes with a small wavenumber k, as 
is usually done when obtaining cosmological constraints, 
both the intrinsic galaxy fluctuations and the lensing fluc- 
tuations are in the linear regime. In real space, even if 
the pair separation is large, one cannot help but mix up 
linear intrinsic galaxy fluctuations with nonlinear lensing 
fluctuations, as long as one considers the LOS orienta- 
tion. Incidentally, this implies that a linear galaxy bias 
is a better approximation in Fourier space than in real 
space. 

The above findings suggest that in precision measure- 
ments of the galaxy power spectrum, such as those that 
attempt to use the baryon oscillation scale to constrain 
dark energy the simplest way 

to immunize against magnification bias could be to go to 
Fourier space, and remove from consideration the small 
fc|l (transverse) modes where magnification bias has the 
largest effect. However, in a photometric redshift survey, 
these are probably the modes with the highest signal- 
to-noise, and the Fourier space also suffers from possible 
complications due to a non-Poissonian shot noise [3J] . A 
better strategy is perhaps to face the magnification bias 
corrections head on, and use the full 3D information to 
constrain and measure them - after all, they contain in- 
teresting cosmological information too. 

Precisely how magnification bias corrections might 
shift the baryon oscillation scale is a subject worthy of 
a separate paper. The precise shift will depend on ex- 
actly how the oscillation scale is extracted from data. A 
preliminary investigation in real space, where the acous- 
tic oscillations manifest as a single local maximum whose 
position is easily defined, was presented in paper I: as dis- 
cussed at the end of miBl the impact on measurements 
of the dark energy equation of state can be significant 
(shifting it by up to ~ 15%). It should also be empha- 
sized that the baryon oscillations are not the only large 



scale features of interest. The radiation-matter equality 
peak at around k ~ 0.01 h/Mpc contains valuable cos- 
mological information, but it can be seen from Fig. [5]-[S] 
panels (b) that the power spectrum at this scale is likely 
significantly affected by magnification bias. 

We have incorporated the effects of peculiar mo- 
tion (and redshift inaccuracy) on the power spectrum 
anisotropy in mil] (the Alcock-Paczynski effect is further 
incorporated in Appendix [X| . The main conclusion is 
that the lensing induced features remain rather robust, 
thanks to their localized nature to small fc|| 's. It should be 
in principle possible to separately measure from data the 
galaxy-galaxy, galaxy-magnification and magnification- 
magnification power spectra, exploiting their different 
shapes in the fc|| — k± plane (see eq. [35]). Exactly 
how to do so in an optimal fashion when faced with 
noisy data deserves further study. In particular, differ- 
ent galaxy types have a different galaxy bias and number 
count slope; how to best weigh their relative contribu- 
tions to the observed power spectrum? 

Lastly, the findings in this paper and paper I cast a 
new light on the well known excess correlations seen in 
pencil beam surveys [H, H^. Could these be the re- 
sult of enhanced correlations due to magnification bias, 
particularly if baryon oscillations are taken into account 
(see [131)? Could the peculiar features seen in the power 
spectrum analysis [11] be due in part to the correspond- 
ing multiplicative windows G and |Vt^||P? We hope to 
address these questions in the future. 
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APPENDIX A: INCORPORATING THE 
ALCOCK-PACZYNSKI EFFECT 

In this Appendix, we incorporate the Alcock-Paczynski 
distortion of the correlation function and power spectrum 
32, 33!| . The observed redshifts and angles are converted 
to comoving radial and transverse distances using the 
Hubble parameter and angular diameter distance. Let us 
denote the parameters used in such a conversion by H^^ 
and x^^- Suppose the true values for these parameters 
are H and x- Then, the two-point correlation function, 
taking into account possible Alcock-Paczynski distortion. 
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is related to the true two-point correlation by 

eb's(xi^P;x2^P) = Cobs(xi;X2) (Al) 

where the right hand side is the true two-point correlation 
given in eq. and Sx^^ = SxH/H^^ and 6xf = 

Ignoring for the moment redshift distortion due to pe- 
culiar motions, putting eq. (fTO|) into eq. (jAip . we have 



(A2) 



^,,{\/6x^^\l-eHr + 6xr\l-e^r) 

+ (1 - e^^)<5x^P/((l - e^)Sxi^) + ^((1 - e^)Sxi^) 



where we have defined 



xlx 



AP 



(A3) 



The important point is this: the exercise of separating 
the three contributions to the observed galaxy correla- 
tion, outhned after eq. (fTU|) (see also Fig. 2 in paper I), 
still works, with minor modifications. For a fixed 5x^ ^ 
the magnification bias corrections (1 — e//)5x'^^/ + 5 still 
dominate over ^gg in the limit of a large 5x^^ ■ This al- 
lows one to fit for the slope (1 — e/f)/ and the intercept 
g. The cosmology dependent factor of 1 — e/f can be 
absorbed into the galaxy bias factor that is present in 
/, which describes the galaxy-magnification correlation. 
The extrapolation back to 5x^^ = gives the intercept g 
which describes the magnification-magnification correla- 
tion. Both / and g are determined up to an overall rescal- 
ing of their argument 5xa_ , which is an uncertainty that 
is always present to the extent the cosmology dependent 
angular diameter distance is uncertain, ^gg can be ob- 
tained by subtracting the inferred galaxy-magnification 
and magnification-magnification contributions from ^^^^ . 

It is common to consider the monopole (eq. 16 in paper 
I). Assuming ejj, <C 1, it can be shown that 



monopole of ^(^^(^a;"^^) = 
igg{Sx^''{l~[eH + 2e^]/3)) 



(A4) 



+ {1 - eff)<5x^^/(<5x^^(l - e^)) + .g(5x^^(l - e^)) 
where 

f{6x^^)= f{6x^^ fim9^)cos9^sm9^d9^ (A5) 

Jo 

g{5x^'^)= / g{5x^^ s\i\ex)s\i\9xd9x 
Jo 

and Ox is the angle between the separation vector and 
the LOS. 

The appearance of the factor 1 — [en + 2ex]/3 in the 
argument of ^gg is the origin of the common statement 
that the baryon oscillation scale measures the combi- 
nation ix^/Hy/^ [14], if one examines the monopole. 



Note, however, the presence of the anisotropic correc- 
tions introduced by magnification bias implies the ob- 
served monopole is no longer related to the true monopole 
by this overall rescaling of Sx. 

Let us give the expression for the observed correlation 
function in the presence of both peculiar motions and the 
Alcock-Paczynski effect: 



(Pk 



(A6) 



b 



e II "e 



/ (xi -x'^)^ (X2 -Xa)^ 



27rcr2 " 

w rdx' 

Jo 



(x~x')x' 



(1 + ^') 



/\2 



-P (7' k , \p^'^^ 



where fo and q are as defined in eq. ([57)1 and 
and Jx"^^ and (5x are related as in eq. (jAip . The first 
term on the right accounts for both the Kaiser effect 
and, in a crude form, the finger-of-god effect due to ei- 
ther virialized motions or photometric redshifts. The sec- 
ond term on the right represents the contribution from 
galaxy-magnification correlation. The finger-of-god dis- 
persion has an effect here, but its effect is negligible for 
large LOS separations (|xi~X2| 3> <7z) where the galaxy- 
magnification correlation has the largest effect (i.e. the 
integral over x'l and x'2 yields roughly |xi - X2I = \5x\)- 
The third term on the right is exactly ^^i^ in eq. ([S]): 
the finger-of-god dispersion has no effect on this term, 
because it has no intrinsic dependence on the LOS sep- 
aration 5x- In other words, to good approximation 
(unless az is very large), the galaxy-magnification and 
magnification-magnification terms take the same form as 
in eq. (|A2[) . This means the exercise of separating the 
three different contributions to the observed correlation 
can be repeated here as well. 

Lastly, let us give the corresponding expression in 
Fourier space. Defining k^^, which is related to k by 



fc||iJ^P/i? 



AP 



kj.x/x 



AP 



(A7) 



the Alcock-Paczynski distorted power spectrum is 
related to the true (windowed) power spectrum Pobs by 



^AP/, apn _ H {x^^ f 



where Pobs is given by 
Pobs(k) = 
d^k' 



HAP ^2 



Pobs(k) 



(A8) 



(A9) 



(27r)3^^^*-^'^''' I ^ 



b 
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+ qb{l + z)GAk\\) 



d^k'j 



Jo 



(x - x')x' 



P,nm{z',k'^)\W^{k^-k'^X'/xW 



It is instructive to integrate out the convolving win- 
dows following eq. and obtain: 



P^^^A^^n = {l + eH){l + e^rP,,{-z,k) 



(AlO) 




q(l + z)G(fc||) 



Jo J^^mmyZ, K) 



where we have adopted the approximation ~ G as 
is done in eq. ([55]) . Note how the distinctive mul- 
tiplicative windows G and |W^||(fc||)P, which are sig- 
natures of the galaxy-magnification and magnification- 
magnification correlations, remain intact in the presence 
of the Alcock-Paczynski effect. They give rise to an en- 
hancement in the observed correlation that is localized 
to low fc||'s. Just as in the case of the real space corre- 
lation function, quantities such as Pgm and Pmm (which 
are related to Fourier transforms of / and g in eq. \A2\ ) 
can be determined from data only up to a rescaling of 
their argument k± (eq. [A7| ) (and in the case of Pgm, 
an overall normalization which can be absorbed into the 
galaxy bias). 
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